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ABSTRACT 

The central part of a dark matter halo reacts to the presence and evolution of a bar. Not 
only does the halo absorb angular momentum from the disk, it can also be compressed and have 
its shape modified. We study these issues in a series of cosmologically motivated, highly resolved 
N-body simulations of barred galaxies run under different initial conditions. In all models we 
find that the inner halo's central density increases. We model this density increase using the 
standard adiabatic approximation and the modified formula by Gnedin et al. and find that halo 
mass profiles are better reproduced by this latter. In models with a strong bar, the dark matter 
in the central region forms a bar-like structure ( "dark matter bar" ) , which rotates together with 
the normal bar formed by the stellar component ("stellar bar"). The minor-to-major axial ratio 
of a halo bar changes with radius with a typical value 0.7 in the central disk region. DM bar 
amplitude is mostly a function of the stellar bar strength. Models in which the bar amplitude 
increases or stays roughly constant with time, initially large (40%-60%) misalignment between 
the halo and disk bars quickly decreases with time as the bar grows. The halo bar is nearly 
aligned with the stellar bar (~ 10° lag for the halo) after ~ 2 Gyr. The torque, which the 
halo bar exerts on the stellar bar, can serve as a mechanism to regulate the angular momentum 
transfer from the disk to the halo. 



Subject headings: galaxies:kinematics and dynamics 
methods:iV— body simulations 



galaxies : evolution 



galaxies: halos 



1. Introduction 

The cold dark matter scenario for galaxy for- 
mation predicts that galaxies should be embed- 
ded within massive, extended, hot, and cuspy dark 
matter halos. The halos can have as much as 
twenty times more mass than the disk. They ex- 
tend well beyond the visible galaxy. Halos are 
hot because they are supported against gravity by 



a large velocity dispersion. Halos are cuspy be- 
cause their central density profiles go as r~ a , with 
a ~ 1.0. Unfortunately, for historical reasons 1 , 
for simplicity, and because of lack of computer re- 
sources, most studies of stellar bars have assumed 
unrealistic halos; from rigid ones (papers in the 



*A well established galactic- scale cold dark matter cos- 
mogony did not arise but until late 1990s. 
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1970s, with few exceptions) to live but small ones 
(papers in the 1980s and 1990s). The need for a 
realistic halo component is being recognized and 
nowadays models in which this dark component 
satisfies some or all of the above requirements are 
more common to find in the literature (Debattista 
& Sellwood 2000; Athanassoula 2003; Valenzuela 
& Klypin 2003; O'Neill & Dubinski 2003; Holley- 
Bockelmann et al. 2005). 

A systematic study of the structural and dy- 
namical changes that the dark matter (DM) com- 
ponent experiences in the presence and during 
the evolution of a bar is missing. Most studies 
that address the issue have focused on the angu- 
lar momentum transfer from the bar to the halo 
(e.g., Weinberg 1985; Combes et al. 1990; Lit- 
tle & Carlberg 1991; Hernquist & Weinberg 1992; 
Athanassoula 1996; Debattista & Sellwood 2000). 
For example, Weinberg (1985) and Hernquist & 
Weinberg (1992), using rigid bars, found that a 
bar loses its angular momentum in few rotation 
periods because of the dynamical friction that the 
dark matter exerts on it. This angular momen- 
tum transfer flattens out the initially cuspy halo 
density profile (see also Holley-Bockelmann et al. 
2005). 

There is no doubt that the disk loses angular 
momentum due to this transfer mechanism, but 
the amount and the rate at which this happens is 
still a matter of debate. Fully self-consistent N- 
body simulations with live disk and dark matter 
halos (e.g., Athanassoula 1996; Debattista & Sell- 
wood 1998,2000; Athanassoula & Misiriotis 2002; 
Valenzuela & Klypin 2003) show that bars slow 
down far less than predicted by Weinberg (1985) 
and Hernquist & Weinberg (1992). Debattista & 
Sellwood (1998, 2000) find in their massive halo 
models, i.e., those for which the contributions of 
the disk and the halo to the mass in the central 
region are comparable, that the disk loses about 
40% of its initial angular momentum in ~ 10 Gyr. 
During this time interval the bar pattern speed, 
fl p , decreases by a factor of four. Valenzuela & 
Klypin (2003), in simulations with much better 
force resolution and a more realistic cosmological 
halo setup, find a decrease in fi p of only a factor 
of two in <~ 7 Gyr (see their model A2). 

The response of the DM halo to the transfer 
of the stellar material from the middle parts of 
the disk closer to the center has been studied re- 



cently by Holley-Bockelmann et al. (2005) and 
Sellwood (2003) in relation to the cuspy problem. 
Cosmology predicts cuspy inner density profiles 
that appear to disagree with those derived from 
rotation curves of dwarf and low surface bright- 
ness galaxies (e.g., de Blok, Bosma, & McGaugh 
2003; Weldrake, de Blok, & Walter 2003), but see 
also Swaters et al. (2003); Rhee et al. (2004). 
Holley-Bockelmann et al. (2005) find a decrease 
in the halo central density, while simulations by 
Sellwood (2003) show the contrary. Here we dis- 
cuss halo density profiles within the context of adi- 
abatic approximation models (e.g., Blumenthal et 
al. 1986; Ryden & Gunn 1987; Gncdin et al. 
2004). Do dark matter halos respond to the in- 
crease of stellar mass in the center by contracting? 
If so, can this increase in density be modeled by 
the standard adiabatic approximation? In agree- 
ment with Sellwood (2003), we will show that 
in our models the halo central density increases, 
sometimes by as much as a factor of 3. 

The influence of the bar on the morphology of 
a halo or a spheroid is a point that has only been 
touched briefly in the literature. Yet, it is an issue 
that deserves a more detailed analysis. It may be 
important for accurate testing of cosmological pre- 
dictions, which are becoming a reality (e.g., Ibata 
et al. 2001). For example, in our bar models 
we find that the initially spherical halo flattens to 
values comparable to those found in cosmological 
halos, for which the minor-to-major axial ratio s is 
s « 0.7 (e.g., Jing & Suto 2002). This correlation 
between the morphology of the stellar component 
and that of the DM halo has been also seen in 
cosmological gasdynamic simulations of galaxies 
(Kazantzidis et al. 2004). 

Returning to bar simulations, Combes et al. 
(1990) found no significant flattening of the bulge 
particles (s « 1) in their experiment A, for which 
Mb/Md = 0.5, where Mb and Md are the mass of 
the bulge and the disk, respectively. Debattista & 
Sellwood (2000) studied the halo response to the 
bar by measuring the second harmonic of the an- 
gular DM particle distribution and found that the 
phase difference between the disk and the halo bar 
gradually goes to zero. This effect was also seen 
in the simulations by O'Neill & Dubinski (2003): 
after t <~ 10 Gyr the orientations of bars (halo 
and disk) stays within 10° of one another. The 
halo bar found by O'Neill & Dubinski (2003) was 
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not very elongated. The axial ratio was only 0.88. 
Halo bars have also seen by Holley-Bockelmann et 
al. (2005, see their Fig. 2) and Athanassoula 
(2005a, b). Based on preliminary results of simu- 
lations with strong bars, Athanassoula (2005a, b) 
finds that halo bars are prolate-like, with inner ax- 
ial ratios being 0.7 or 0.8. Contrary to results of 
other groups, her halo bars are roughly aligned to 
the disk ones at all times. 

More recently, in order to understand the role 
of the halo in the formation and evolution of the 
bar, Athanassoula (2002, 2003) has carried out 
an orbital-resonant study of the disk and halo 
particles in models in which the disk parameters 
are kept fixed while varying the degree of the 
halo influence. She found significant differences 
in the halo orbital properties between her halo 
and disk-dominated models. The inner Linblad 
resonance, ILR, is present in the halo-dominated 
models, while in the disk-dominated ones is not. 
In general, she found that there is much more ma- 
terial between the ILR and corotation in the for- 
mer models than in the latter ones. Because these 
halo-dominated models have stronger bars than 
their counterparts- disk-dominated ones, this re- 
sult shows the positive effect that a live halo has 
on the bar growth. 

The paper is organized as follows: In section 
2 we introduce the models to be explored in this 
paper. In Section 3 we present the evolution of 
the amplitude and pattern speed of the bar. Sec- 
tion 4 is devoted to an analysis of models of adi- 
abatic compression, which are used in the context 
of the inner disk density growth. The influence of 
the disk and stellar bar on the initially spherically 
symmetric halo is discussed in section 5. In partic- 
ular, we will show that once formed, the halo bar 
tends to rapidly align with the stellar bar. A brief 
discussion of the redistribution and evolution of 
the disk angular momentum is given in section 6. 
In section 7 we discuss our results and summarize 
our main conclusions. 

2. Models and simulations 
2.1. Initial conditions 

The initial conditions are described in Valen- 
zuela & Klypin (2003). Here we briefly summa- 
rize parameters of our models. The system of a 
halo and a disk, with no initial bulge or a bar, is 
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Fig. 1. — Initial circular velocities profiles for 
models D cs and K c b (top panels), and D c \, and 
A w b (bottom panels). Total, stellar, and halo com- 
ponents are denoted by solid, dotted, and dashed 
lines, respectively. In the disk-dominated models 
(top panels) the system evolves secularly produc- 
ing a weak bar at later stages, while in the halo- 
dominated models (bottom panels) a strong bar is 
present after 5 Gyr of evolution. 



generated using the method of Hernquist (1993). 
In cylindrical coordinates the density of the stellar 
disk is approximated by the following expresion: 



Pd(R, z) 



M d 



e- R/Rd sech 2 (z/z Q ), (1) 



where Md is the mass of the disk, Rd is the scale 
length, and zq is the scale height. The latter is 
assumed constant through the disk. The radial 
and vertical velocity dispersion are given by 



a(R) = Q 



3.36GE(i?) 
k{R) 



a 2 z (R) = ttGzq^(R), 



(2) 

where k is the epicyclic frequency, and Q the sta- 
bility or Toomre parameter. Our models keep Q 
fixed along the disk. The azimuthal velocity and 
its dispersion are found using the asymmetric drift 
and the epicycle approximations. 

The models assume a NFW density profile 
(Navarro et al. 1997) for the halo component, 
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Table 1 
Simulation Parameters 



Parameter 


A w b 


D cs 




D hs 


K cb 


K w b 




Disk mass (10 10 M Q ) 


4.28 


5.0 


5.0 


5.0 


5.0 


5.0 


5.0 


Total mass (10 12 M ) 


2.04 


1.43 


1.43 


1.43 


1.43 


1.43 


1.43 


Disk exponential length i?d(kpc) 


3.50 


2.57 


3.86 


2.57 


3.86 


3.86 


3.86 


Disk exponential height zo(kpc) 


0.14 


0.20 


0.20 


0.20 


0.20 


0.20 


0.20 


Toomre parameter Q 


1.6 


1.3 


1.3 


1.8 


1.2 


1.5 


1.8 


Halo concentration C 


15 


17 


17 


17 


10 


10 


10 


Number of disk particles (10 5 ) 


2.0 


2.33 


2.33 


4.60 


2.33 


2.33 


2.33 


Total Number of particles (10 6 ) 


2.7 


2.5 


1.9 


3.80 


2.3 


2.7 


2.3 


Particle mass (10 5 M ) 
Time-step (10 4 yr) 


2.14 


1.07 


2.14 


1.07 


2.14 


2.14 


2.14 


1.5 


2.0 


1.5 


1.5 


1.2 


1.5 


2.0 


Minimum cell size (pc) 


22 


48 


22 


22 


24 


22 


48 


Disk-to-DM mass ratio inside Rd 


0.82 


2.03 


0.98 


2.03 


1.95 


1.95 


1.95 


Strong Bar after 5 Gyr(A 2 > 0.4) 


yes 


yes 


yes 


no 


no 


no 


no 



which is described by 
Pdm{t) = 

„3 



Po 



M mr = inpori 



ln(l + c) 



1 + c 



x = r/r s , (3) 

Rvir / ,\ 

c =—> ( 4 ) 



where M V i r , R V i r , and c are the virial mass, the 
virial radius, and concentration of the halo, respec- 
tively. Given M vir , the virial radius is found once 
a cosmology is adopted 2 . The equation (4-56) of 
Binncy & Trcmaine (1987) and the assumption 
of isotropy in the velocities allow us to determine 
the radial velocity dispersion as 



' r.DM 



PDM J r 



GM(r) , 
Pdm 5 — dr, 



(5) 



where M(r) is the mass contained within radius r 
and G the gravitational constant. 

The dark matter halo is truncated at the virial 
radius, which for our D and K models is 287 kpc, 
while for our A w b model it is 323 kpc (see Table 
1). The disk component is realized with particles 
of equal mass. The halo is composed of particles 
of different mass placed inside out as we go from 
less to more massive particles. The smallest par- 
ticles have a mass equal to those of the disk. The 



2 We adopt the fiat cosmological model with a non- vanishing 
cosmological constant with Qq = 0.3 and h = 0.7. 



multiple-mass scheme is designed to reduce the to- 
tal number of particles and, thus, the CPU time. 

Parameters of our models are presented in Ta- 
ble 1. They were chosen to cover a relatively 
wide range of the parameter space: there are hot 
and cold models, Khb versus K c b] there is a disk- 
dominated D cs and a sub-maximum D c b model. 
There are models with different DM halo con- 
centrations: D c b and K c t,. The parameter values 
are similar to those used by Valenzuela & Klypin 
(2003), whose models were chosen to mimic the 
properties of the Milky Way galaxy (Klypin ct 
al. 2002). More specifically, models labeled with 
D have all the same disk mass, total mass, and 
halo concentration. This latter is above the me- 
dian ~ 13 for a halo of this mass and cosmology 
adopted here (Bullock et al. 2001a). On the other 
hand, models K differ one from other only by the 
value of the disk stability parameter Q. 

As a measure of how dense is the halo in the 
central region we show in the Table 1 the disk-to- 
DM mass ratio evaluated at the disk exponential 
length Rd- The values are about 1-2. In other 
words, from the beginning the dark matter is sub- 
dominant in the central region. This is broadly 
consistent with the mass modeling of the Milky 
Way galaxy (Klypin ct al. 2002): even for cosmo- 
logically motivated cuspy DM halos, the density 
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Fig. 2.— Evolution of the bar pattern speed (first and third rows) and the bar amplitude (second and fourth 
rows) for six of our seven models. Models D CSl K w b, and K^b evolve very slowly after the bar formation. 
In models Dh s and K c b (not shown in the plot), the bar amplitude, after reaching a peak, decreases as a 
function of time. In models A w b and D c b, at late stages of the evolution, the bar develops a second growth 
period. The dashed curves in the Khb model show results obtained with the GADGET code. 
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of the dark matter must be smaller than the den- 
sity of the baryons in the central region of galax- 
ies. This is different from what is often assumed 
in simulations of barred galaxies. For example, in 
DM dominated models of Athanassoula & Misiri- 
otis (2002) the halo mass dominates everywhere 
inside 2Rd- This is not consistent with cosmologi- 
cal predictions. We do allow variations of the DM 
mass in the central region. For example, the model 
D c b has a denser halo as compared with models K . 
Models D c b and A w b, in which initially the disk is 
less dominant, develop strong bars. The last row 
of Table 1 classified the models according to stel- 
lar bar strength after 5 Gyr of evolution. We label 
them as strong bar models if they develop a bar 
with an amplitude larger than 0.4. 

Figure 1 shows the initial circular velocities pro- 
files V c = y/GM(r)/r for all the components. 
This estimate of V c assumes spherical symme- 
try. The circular velocity estimated in the plane 
of the disk using the real non-spherical distribu- 
tion V c ,z=o = \J g(r)r, where g(r) is the gravita- 
tional acceleration in the disk plane, is very close 
(whithin 5 percent) to the circular velocity V c at 
radii larger than 1 kpc. At smaller distances it is 
slightly (10-20%) below V c . 

2.2. Numerical simulations 

The simulations were run with the Adaptive 
Refinement Tree (ART) code (Kravtsov et al. 
1997). The ART code starts with a uniform grid, 
which covers the whole computational box. This 
grid defines the lowest (zeroth) level of resolution 
of the simulation. Two grids with 256 3 and 128 3 
cells in the zeroth- level and boxes of 1.43 Mpc and 
1.0 Mpc across, respectively, were used by the sim- 
ulations presented here. The models are placed at 
the center of the grid and far from the periodic 
images. Because the size of the models are small 
as compared with the box sizes, the effects of the 
periodical images are negligible. For example, at 
a distance of 100 kpc from the center of a model, 
in the 1.43 Mpc box, the relative contribution of 
periodic images is less that 7 x 10~ 7 of the main 
galaxy force. Even at 250 kpc, which is close to 
the virial radius for the models, the tidal force is 
only 6 x 10 -5 of the force from the central im- 
age. Bar formation and evolution are practically 
unaffected by the periodic boundary conditions. 

The ART code achieves high spatial resolu- 



tion by refining the base uniform grid in all high- 
density regions with an automated refinement al- 
gorithm. The standard Cloud-In-Cell algorithm 
is used to compute the density and gravitational 
potential on the zeroth-level mesh with periodi- 
cal boundary conditions. The code then reaches 
high force resolution by refining all high density 
regions using an automated refinement algorithm. 
The refinements are recursive. A refined region 
can also be refined. Each subsequent refinement 
level has half of the previous level's cell size. This 
creates a hierarchy of refinement meshes of differ- 
ent resolution, size, and geometry covering regions 
of interest. Because each individual cubic cell can 
be refined, the shape of the refinement mesh can 
be arbitrary and effectively match the geometry 
of the region of interest. This algorithm is well 
suited for simulations of a selected region within 
a large computational box, as in the simulations 
presented below. 

The criterion for refinement is the local den- 
sity of particles. If the number-density of parti- 
cles in a mesh cell (as estimated by the Cloud- 
In-Cell method) exceeds the level nthresh, the cell 
is split ("refined") into 8 cells of the next refine- 
ment level. The refinement threshold depends on 
the refinement level. The threshold for cell refine- 
ment was low on the zeroth level: n t hrcsh(0) = 2. 
Thus, every zeroth-level cell containing two or 
more particles was refined. The threshold was 
higher on deeper levels of refinement: n t hresh = 3 
and nthresh = 4 for the first level and higher levels, 
respectively. Note that the code actually does not 
count the number of particles in a cell to decide 
whether the cell should be split or not. It uses 
the density field, which is less noisy than count- 
ing particles. It then filters the density field to 
make the map of refinement even smoother. On 
average the algorithm maintains 1-4 particles per 
cell. Yet, some cells have few particles in them and 
some can be even empty. Those cells are not used 
for estimates of the gravitational acceleration act- 
ing on particles. For each particle the code checks 
whether all 56 cells surrounding the cell hosting 
the particle are on the same level of refinement. 
If they are not, then the code assigns the particle 
to a lower resolution and checks again for the re- 
finement level of all its neighbors. In two of our 
models (D cs and K^b) we limited the resolution to 
not more than 7 levels of refinement. 
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The resolution changes from place to place and 
it also changes with time. During the initial stage 
of evolution - the formation of the bar - the den- 
sity in the central part increases substantially. For 
example, in the model K c b the density at 250 pc 
increases by a factor of 10 by the end of the evo- 
lution. The code reacts to this density increase 
by changing the number of high-resolution cells. 
The lion's share of this increase happens during 
the first 0.5 Gyr when the system develops very 
strong non-axisymmetric fluctuations. For exam- 
ple, in the case of the model K c b the density at 
the central 250 pc increases 6 times during first 
t = 0.5 Gyr. It will take another 5 Gyr to double 
it. The code increases the number of high resolu- 
tion 24 pc cells 14 times during the first 0.5 Gyr of 
evolution. The number of cells then changes only 
by 30 percent during the rest 5 Gyrs. The increase 
on lower 48 pc level of resolution is even smaller: 
the number of cells increased only by 20 percent 
during the whole evolution. Note that this implies 
that the size of the high resolution region changes 
only by 6 percent. In other words, the increase of 
the resolution with time is very mild and most of 
it happens before the bar forms and settles. 

The spatial resolution varies with the distance. 
Typically in the models the central roughly ellip- 
tical region with large axis 1.5 kpc (4 kpc) is re- 
solved with the 20-25 pc (40-50 pc). The whole 
disk (up to radius of about 15 kpc) is resolved 
with 80-100 pc cells. The force resolution in our 
simulations is twice larger than a cell size. In- 
side one resolution element (roughly a sphere of 
two cell radii) the code has on average 64 parti- 
cles. This large number of particles implies that 
the code does not produce too high resolution 
for the number of particles present in our mod- 
els. The effect of close encounters is actually even 
smaller than what it seems from naive counting 
particles. The natural scale of the close encoun- 
ters is Chandrasekhar's minimum impact param- 
eter 6 m ; n = Gm/V 2 , which is the distance for 90 
degree deflection (here m is the particle mass and 
V is the relative velosity). One does not want to 
get close to it, but it still gives a scale for the 
small-distance scattering. In the central 1 kpc re- 
gion of our models the rms velocity of particles 
is about 200 km/s. With the mass of a particle 
of 2 x 10 5 M Q , we get 6 min = 2 x 10~ 2 pc. This 
is 2000 times smaller than our resolution of 50pc! 



And this still overestimates the effect because we 
have many particles inside the resolution element. 
At larger distance, say 5 kpc, the rms velocity in 
the disk is smaller. Again, taking a typical rms 
velocity from our models, we get 6 m i n = 0.1 pc. 
With the cell size of 88 pc (resolution 160pc) we 
are still 2000 times away from 6 m i n . To summa- 
rize, effects of close collisions are definitely small 
in our models. 

During the integration, spatial refinement is ac- 
companied by temporal refinement. Namely, par- 
ticles moving on each level of refinement are ad- 
vanced using the leap-frog scheme. When a par- 
ticle moves to another level, its velocity is re- 
interpolated. The time step decreases by factor 
two with each refinement level. This variable time 
stepping is very important for accuracy of the re- 
sults. As the force resolution increases, more steps 
are needed to integrate the trajectories accurately. 

2.3. Code tests 

ART code has been extensively tested during 
the last ~ 10 years. A list of tests and details are 
presented in Valenzuela & Klypin (2003). Here 
we give a brief summary of those which are im- 
portant for our problems. Among other tests, 
Kravtsov et al. (1997) tested the code using the 
spherical accretion model. A small initial seed is 
placed in a homogeneous expanding background. 
Shells surrounding the seed start to accrete on 
the seed, which results in the built up of a very 
centrally concentrated object. An analytical so- 
lution for this complicated process is known (the 
Bertschinger solution). In the test, the density in- 
creases by almost 5 orders of magnitude. At the 
end, no detectable deviations of the numerical so- 
lution provided by ART and the analytical solu- 
tion were found. This is an important test because 
it shows that the code accurately treats collapsing 
systems. In the simulations of barred galaxies, this 
collapse happens during the initial stage of forma- 
tion of bars (first ~ 1 Gyr of evolution) , when the 
density in the central region increases few times. 
Asymmetric collapse was tested (Kravtsov et al. 
1997) using the Zeldovich approximation (plane 
wave collapse). 

Valenzuela & Klypin (2003) present additional 
tests, which focus on a long-term stability of equi- 
librium systems. These tests are relevant for later 
stages of evolution of barred models, when the sys- 
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tern changes its structure very gradually. Valen- 
zuela & Klypin (2003) used 200,000 particles 
to set equilibrium for a high concentration NFW 
halo. The test is actually very difficult for the 
code. Orbits of particles in the NFW halo are 
very elongated. Thus, particles move again and 
again from high density regions with high force 
resolution to outer regions with low resolution and 
back. Any numerical defects at boundaries of re- 
gions of different resolution would result in evolu- 
tion of the system or in excessive scattering of tra- 
jectories. No defects were detected for gigayears 
of evolution. Using trajectories of individual par- 
ticles, Valenzuela & Klypin (2003) estimated ef- 
fects of the two-body scattering in the simulations. 
The two-body scattering is clearly present. It is 
consistent with predictions of the Chandrasekhar 
approximation (after some scaling due to inhomo- 
geneous distribution of density). Yet, the time for 
the two-body scattering is very long: for a typical 
run with 3 x 10 6 particles the two-body relaxation 
time is rts 4 x 10 4 Gyr. So, it is practically negli- 
gible. 

Numerical two-body scattering is often dis- 
cussed in conjunction with effects of resonances 
in barred galaxies (e.g., Weinberg, & Katz 2005). 
Excessive scattering may prevent accurate treat- 
ment of resonances, which are important factors 
in disk galaxies. Formation of bars in dark matter, 
discussed below, is an example of resonant inter- 
action between stellar bars and the dark matter. 
Ceverino & Klypin (2005) studied resonances in 
some of our models. A large number of strong 
resonances were detected, including the familiar 
co-rotation, the inner and the outer Linblad res- 
onances. Thus, the two-body scattering in our 
simulations is small enough for the resonance to 
be present and to be very important factors in the 
evolution of the systems. 

After reaching a peak in their amplitude, bars 
in our models can either grow, gradually decline, 
or remain in a steady state with little change in the 
amplitude. For models with relatively low central 
densities (e.g., A w \> and Khb) the evolution with 
time is robust against changes of the numerical 
parameters: the timestep, the number of the par- 
ticles, and the force resolution. We rerun some 
of our models doubling the number of particles 
while keeping the same small timestep and the re- 
sults do not change significantly. We also run the 



K hb model with the GADGET code (Springel et 
al. 2001). In this case we used only 10 5 parti- 
cles in the disk and 10 6 total. The force resolu- 
tion was limited to 100 pc (Plummer softening). 
At the beginning of the simulation the time-step 
was <~ 1.5 x 10 5 yr while by the end of the evo- 
lution it had reduced to ~ 4.4 x 10 4 yr. Figure 2 
shows the evolution of the pattern speed Q p and 
the bar amplitude A 2 for the model Khb run with 
the GADGET and the ART codes. Overall the 
agreement with the ART code is reasonably good. 
After 5 Gyrs of evolution the bar amplitudes and 
the bar lengths are within 10 percent. The pattern 
speeds deviate not more than 20 percent. Note 
that both simulations of the Khb model show the 
same qualitative behavior: the pattern speed is 
barely changing over 5 Gyrs of evolution and the 
amplitude of the bar is gradually decreasing with 
time. 

3. The bar pattern speed 

Figure 2 shows the evolution of the bar pat- 
tern speed ftp and the amplitude of the bar A 2 
for six models. To estimate ft p we first determine 
the orientation of the bar computed applying it- 
eratively the method of the tensor of inertia in 
the plane of the disk (see section 5 below). Q p 
is obtained subsequently by numerical differentia- 
tion: Q p = d(f>/dt, where <p is the position angle 
of the bar. In practice, we use about ten consec- 
utive snapshots for which the increasing function 
<j) is available and make a least-squares fit. Then 
Vl p is given by the slope of the straight line. The 
bar amplitude A 2 is computed as in Valenzuela & 
Klypin (2003). The curves are smoothed out in 
time by using a top-hat kernel. Models show what 
seems to be a generic feature of bar simulations: 
an increase in A 2 is accompanied with a decline 
in Q p and vice versa. A steady amplitude, on the 
other hand, produces a roughly constant £l p . 

4. Adiabatic contraction 

The structure and the dynamics of the halo is 
altered by the presence and evolution of the disk 
and vice versa. For example, halo's inner shape 
changes from the assumed initial spherical con- 
figuration to a triaxial one. We will see that in 
models, in which a strong bar developes, the halo 
becomes prolate. The m — 2 structure of the halo, 
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Fig. 3. — The cumulative halo mass profiles for 
models D cs and D c b. Initial profiles are shown by 
dotted curves. The profile at the last simulated 
epoch t ~ 6 Gyr is presented by full curves. Re- 
sults of the standard and the modified adiabatic 
approximations are presented by the short-dashed 
and long-dashed curves. In the small panels we 
also show the relative mass differences between the 
adiabatic approximations and the mass profiles as 
measured in simulations. 

the DM bar, which arises because of the bar insta- 
bility of the disk, couples with the stellar bar. It is 
an important contributer to the stellar bar brak- 
ing (see section 6 below) . In this section, we study 
the changes in its mass structure. 

It is known (e.g., Athanassoula & Misiriotis 
2002; Valenzuela & Klypin 2003; Avila-Reese et 
al. 2005) that as the stellar bar grows, the inner 
disk density increases. Is this density growth ac- 
companied with a corresponding central halo mass 
increase? If so, can it be modeled by the standard 



S 0.1 



L I I | I I I I I 11,1.1. 

• X : 

V ' 

v 

: >>. : 

: d X 

cs X 


i . 1 . i i . 1 , i i | i i i i 
^^^^ ^ 

; ^X ; 

: K cb X, : 

: 1 1 | INI 1 |X : 


r i i I mi i I i 1 1 r 

: \ : 
: \ : 

x ; 

'X. 

"X 

X 

. D . X 

cb X 

? , , , 1 , . , , : 


: K hb X x ^ ; 



0.5 1 5 10 0.5 1 5 10 

r(kpc) r(kpc) 



Fig. 4. — Spherically averaged DM density pro- 
files for models D cb , Khb, D cs and K cb . Curves 
are coded as in previous figure: initial profile (dot- 
ted curve), the profile at the end of the evolution, 
t ~ 6 Gyr, (solid curve), and the adiabatic ap- 
proximation of Gnedin et al. (long-dashed curve) . 
Noise in initial and final density profiles was re- 
duced by averaging five consecutive timesteps. For 
comparison in the top-right panel we show power 
laws with the slope —1 and —2. 

adiabatic contraction approximation? This latter 
model assumes spherical symmetry, homologous 
contraction, particles in circular orbits, and con- 
servation of the angular momentum (Blumenthal 
et al. 1986; Ryden & Gunn 1987; Gnedin et al. 
2004). Under these assumptions, the final halo cu- 
mulative mass profile can be computed given the 
initial DM and barionic mass profiles, -Mdm(?"), 
Mft(r), and the final barionic profile, Mf,(r/): 

[M DM {r) + M b {r)\ r = [M DM (r) + M b {r f )\ r f . 

(6) 

Given rf and, thus, Mj(rj), we solve the above 
equation for r iteratively. The initial mass profiles 
as well as the final baryonic profiles are taken from 
simulations. To reduce noise, both final and initial 
profiles are smoothed out over few close timesteps. 
Gnedin et al. (2004) tested recently the model in 
cosmological simulations and find that it system- 
atically overestimates the central density. Notice 
that since the orbits of DM particles are on average 
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eccentric the assumption of circular orbits of the 
standard formalism breaks. Gnedin et al. (2004) 
show that the model can be improved if the adia- 
batic invariant is changed from M(r)r to M(f)r, 
where r and f are the current and orbit-averaged 
particle positions, respectively. This modification 
approximately accounts for the eccentricity of par- 
ticle orbits. It was found that r can be described 
by a power law: 

x = Ax w , x = r/r vir , (7) 

with small variations in the A and w parameters. 
We will test the new model adopting mean values 
given by Gnedin et al. (2004): A = 0.85 and 
w = 0.8. 
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Fig. 5. — The disk to dark matter ratio for models 
D c b, Khb, D cs , and K c b- The dotted curves repre- 
sent the initial models. The final models are shown 
with solid curves. The increase in the disk central 
density is produced at the expense of the density 
in middle region. This effect is slightly more pro- 
nounced in in the cold- and disk-dominated models 
D cs , and K cb . 

Figure 3 shows the initial and the final halo 
mass profiles for models D cs (top panel) and D c b 
(bottom panel). The mass profile in the model 
D c b barely changes, while in the model D cs the 
mass at 0.5 kpc, for example, has increased sub- 
stantialy - by a factor of 2.5. Overall, the cen- 
tral mass increase seem to be reproduced better 



by the modification introduced by Gnedin et al. 
(2004). The standard formula overestimates the 
central density, which agrees with the results of 
Gnedin et al. for cosmlogical halos. 

Figure 4 shows the dark matter density profiles 
of models D cs and K c b (cold and disk-dominated 
models), and models Khb and D c b (hot and dark- 
dominated models). The adiabatic density profiles 
are computed by differentiating the correspond- 
ing mass profiles. We decided to show only the 
predictions by the adiabatic contraction formal- 
ism developed by Gnedin et al. (2004) because 
the standard adiabatic contraction formula does 
not produce better results. In order to reduce the 
shot noise, we average five consecutive snapshots 
in both final and initial profiles. As was already 
noticed in the cumulative mass profiles (Fig. 3), 
the increase in the inner halo density is higher in 
disk-dominated and cold models (D cs versus D c b 
or K c b versus Khb)- The inner slope is steeper 
than 7 = —1 in all models except in the model 
D c b- The rather flat inner density slope of the 
model D c b (see left-bottom panel) is related with 
the fact that the inner density profile of the stellar 
disk is also flat in this model. 

To summarize, we find that the increase in the 
halo central density is compatible with the adia- 
batic compression. In particular, we find that the 
inner halo density increase is better reproduced by 
the formalism developed by Gnedin et al. (2004) . 
In models in which this density increase is smaller 
(model D c b, for instance), the adiabatic compres- 
sion tends to overestimate the DM density in the 
central region. 

Figure 5 gives additional information on the 
changes in the dark matter. It shows the ratio 
of the disk to the halo densities as a function of 
radius. At the end of the evolution in all mod- 
els the stellar component strongly dominates the 
central ~ 1 kpc region. 

5. Shape of the dark matter halo 

The ellipticity of different components (DM or 
stellar) as well as the directions of the principal 
axes are determined by iteratively diagonalizing 
the inertia tensor. We start by finding the tensor 
of inertia for all particles inside a spherical shell of 
given radius R. We can take either stellar or DM 
particles depending on what component we study. 
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Fig. 6. — Distribution of disk (top panel) and dark 
matter particles (bottom panel) inside a box of 
10 kpc across for model at t ~ 6Gyr. The 
figure shows the face-on projection of the parti- 
cles. The halo bar is clearly seen in the bottom 
panel. To improve the contrast, we have color- 
coded particles on a gray scale according to their 
local 3D density (a pgplot code kindly provided by 
A. Kravtsov) and plotted only DM particles with 
the z- component of velocity \v z \ < 100 km/s. 

We then find principal semi-axes a > b > c, and 
angles of the tensor. In the next iteration, we find 
the modified inertia tensor for particles inside 



an elliptical shell with the orientation given by the 
inertia tensor on the previous iteration and with 
the semi-axes R, (b/a)R, and (c/a)R: 



dl 



2 , ~2,fc ^3,k 
x l,k T o ' 2 ' 

b/a, s = c/a. 



(8) 

(9) 
(10) 



Here the sum is taken over all particles inside the 
elliptical shell. In equations (8-10) Xi^ (i = 1, 2, 
and 3) are the coordinates of the particle k with 
respect to the center of the disk-halo system and 
dk is the elliptical distance; s is the short-to-long 
axis ratio, and q is the intermediate-to-long axis 
ratio. The center of the system is found iteratively 
as the center of mass of a sphere containing max- 
imum mass. The radius of the sphere is equal to 
the disk scale length Rd- The elliptical bins are 
chosen so as they have roughly the same number 
of particles and are less than 1 kpc width. As ex- 
pected, the minor axis lies almost perpendicular 
to the plane of the disk. 

Figure 6 highlights the presence of the stellar 
and the DM bars in the model D c i,, which has 
a strong bar. The stellar bar shown in the top 
panel has a rectangular-like shape, while the halo 
bar shown in bottom panel is more round. Notice 
that both bars have their major axes pointing to 
approximately the same direction. 

Figure 7 shows the flattening, / = 1 — s = 
1 — c/a, and the ellipticity of the dark matter in 
the plane defined by major and intermediate axes, 
££> = 1 — 5=1 — b/a, as a function of distance R 
to the center for models D cs and D c b at the last 
simulated epoch ~ 6 Gyr (hereafter / and ed arc 
referred to as ellipticities). The plot shows that 
DM bars are typically triaxial. In models with a 
strong stellar bar, as in the model D c b, the DM 
bar can be more prolate than oblate. As shown 
in Figure 7, axial ratios increase with increasing 
radius (thus, the system is getting more round). 
For example, at 15 kpc models D cs and D c b show 
comparable small ellipticities, while in the central 
regions the models have very different flattening. 

Evolution of the DM ellipticities with time is 
shown in Figure 8. By comparing Figure 2 with 
Figure 8 we find similarities of behavior of the 
stellar bar amplitude with the evolution of DM 
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R(kpc) R(kpo) 

Fig. 7. — Flattening / = 1 — c/a (solid curve) and 
ellipticity in the plane of the disk en = I — b/a 
(dashed curve) of the DM distribution, where a, b, 
and c are the major, intermadiate, and the minor 
axes. We show models D cs (left panel) and D c b 
(right panel) at the last measured epoch t ~ 6 Gyr. 
The differences between the weak bar model D cs 
and the strong bar model D c b are clearly seen. 

bar ellipticity. Indeed, in models Dhs and Khb 
the maximum in the stellar bar amplitude Ai and 
the maximum of ellipticity eo are reached approx- 
imately at the same time. Moreover, in the model 
D c b the second period of the stellar bar growth is 
mirrowed by an increase in the ellipticity of the 
DM bar. The gradual decline of A 2 in models 
Dh s , K c b, and to smaller degree in the model Khb 
is imitated by a decline in en- 

At the onset of the evolution, the sole presence 
of the disk produces a flattening of the DM halo 
in the central region. The subsequent evolution of 
the shape of the DM halo depends on the model. 
For example, in model Dhs the ellipticities first 
increase as the stellar bar growths. Then they 
decline maintaining a figure of constant triaxiality. 
On the other hand, in the model D c b the flattening 
/ stays approximately constant after the initial 
period of stellar bar growth. This latter effect, 
along with the growth of ellipticity eo, produces 
in models D c b and A w b (not shown in the figure) 
a near prolate central DM halo by the end of the 
evolution. In summary, Figures 2, 6, and 8 show 
a clear indication of the disk-halo coupling (for 
details see section 6). 

5.1. Bar orientation 

We use the position angle, PA, of the stellar 
bar and the orientation of the DM bar measured 
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Fig. 8. — Evolution of the DM ellipticities, / (top) 
and eo (bottom) for four of our seven models. El- 
lipticities are evaluated at r = Rd, where Rd is the 
scale length of the corresponding model. 

at r = Rd to build the evolution of the bars phase 
difference, Ac/). Figure 9 shows the evolution of A(j> 
for four of our seven models. At Rd the bar align- 
ment within ~ 10° is reached soon after the initial 
stellar bar growth (the DM bar trails the stellar 
bar). The start of the second period of stellar 
bar growth in model D c b coincides with a small 
increase in Acj) - the DM bar gets a larger lag. 
This renewed disalignment quickly dies out and 
the bars again rotate in phase. More significantly, 
this second period coincides with the period of 
a strong influence of the stellar bar on the halo, 
which is reflected in the more extended bar aug- 
ment (see top- left panel, dot-dashed curve). These 
two periods of stellar bar growth and bars align- 
ment also coincide with the two periods in which 
a higher angular momentum lost-rate is observed 
(see Figures 11 and 12). 

The phase difference A<fi of the bars is related 
to the torque experienced by the stellar bar: the 
torque slows down the bar pattern speed. In the 
model A w b (see Figure 2) the oscillations seen in 
Q p , and to some degree also in A 2 , are very likely 
coming from the strong oscillations seen in A</> 
in this model. The onset of the late period of 
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bar growth coincides with a decline of about 40° 
in A</>: the m = 2 halo structure is quickly re- 
oriented toward the direction of the stellar bar. 
Moreover, the rapid growth of the stellar bar in 
the model D cs (see Figure 2) produces a DM bar 
that very soon aligns with the stellar bar. 
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Fig. 9. — Evolution of the phase difference be- 
tween the stellar and the DM bars for four mod- 
els at radius Rd (solid curves). For the strong 
bar model D c b we also show the evolution of A</> 
at other radii: R — Q.5Rd (dashed curve), l.hRd 
(dotted line), and 2.0Rd (dot-dashed curve). The 
curves show that the phase difference is an increas- 
ing function of radius. Interestingly, the two peri- 
ods in which A(j> drops (see solid and dot-dashed 
curves) seem to coincide with the two periods in 
which a higher angular momentum lost-rate is ob- 
served. At R = Rd an alignment within 10° is 
observed in models that end up with strong and 
moderate stellar bars by the end of evolution. 



6. Angular momentum 

The bar growth is intimately related to its slow- 
ing down. The rate of change of £l p correlates with 
the initial conditions. For example, the A-i maxi- 
mum is reached earlier in disk-dominated models 
(compare D cs with D c b), or in cold disk models 
(K c b versus Khb)- During the bar braking the an- 
gular momentum is transfered from the disk to the 
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Fig. 10. — Distribution of the z-component of 
the angular momentum of the stellar particles for 
models K c b (top panel) and D c b (bottom panel), 
shown at different epochs: initial (long-dashed 
curve), at t ~ 1.6 Gyr (dotted curve), and at 
t ~ 5 Gyr, (solid curve). Particles with low 
(high) angular momentum are preferentially at 
small (large) radii. The distributions are quali- 
tatively similar to those shown by Valenzuela & 
Klypin (2003) (see their Figure 12): a peak at 
low L z , which corresponds to the bar, a valley at 
intermediate values of L z , and an excess of of par- 
ticles with very large angular momentum. Notice 
that in the weak bar model K c b the distribution 
barely changes for the period 1.6 — 5.0 Gyr, while 
the redistribution of angular momentum continues 
in the strong bar model D c b- 



halo. 

Figure 10 shows the distribution of the disk 
L z component of the angular momentum for two 
models, the strong bar model D c b and the cold 
model K c b, one of the models with a strong decline 
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in bar amplitude. We show the distribution of L z 
at three moments of time: initial, at ~ 1.6 Gyr, 
and at ~ 5 Gyr. Disk particles in the intermediate- 
L z zone lose angular momentum and move to the 
low-L z region close to the center. At the same 
time, the number of particles with very high-L z , 
which are preferentially at larger radii, increases. 
Note that the number of intermediate-L z particles 
in the model K c b barely changes from t ~ 1.6 to 
~ 5 Gyr, while this is not the case in the model 
D c i,. This behavior is compatible with the much 
higher angular momentum loss by the model D c b 
during this period of time shown in Figure 1 1 . 

The evolution of the relative change in L z is 
plotted in Figure 11, where we show four models: 
three cold models D cs , K c b, and D c b and one hot 
and disk-dominated model Khb- In the K c b model, 
the amount of angular momentum lost by the disk 
is very small during the declining phase of the bar 
amplitude. Moreover, the disk in models D cs and 
D c b has lost, by the end of the evolution, about 
10% of its angular momentum. 
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Fig. 11. — Evolution of the disk angular momen- 
tum for four models. We show two models, D c b 
and D cs , that end up with strong bars at the end 
of the simulated period and two other models, K c b 
and Khb- Model D c b presents an inflexion point or 
zone, in which the rate of angular loss changes its 
decreasing trend to an increasing one. This time 
coincides with the second period of bar growth. 



It is interesting to estimate what fraction of the 
torque experienced by the stellar disk is coming 
from the interaction between the dark matter bar 
and the stellar bar. Unfortunately, the torque is 
difficult to estimate accurately because of the com- 
plex structure of the bars: the bars change shape, 
mass, and length as the evolution proceeds. Still, 
we can start with simple scaling relations. The 
torque should be proportional to masses and radii 
of the bars. It is also proportional to the strength 
of the bars, which can be measured by ellipticities 
q = 1 — b/a, where a and b are the major and mi- 
nor semi-axis in the plane of the disk. The angle 
A<f> between the bars also affects the torque. From 
the geometry of the problem, we expect that the 
torque is zero, when the angle is either zero or 90°. 

In order to test the scaling relations, we esti- 
mate the torque between two elliptical-shape ob- 
jects, which roughly mimic our bars. Boundary 
of each object is assumed to be a triaxial ellip- 
soid with semi-axes a > b > c: (x/a) 2 + (y/b) 2 + 
(z/c) 2 = 1. For the ellipsoid representing the stel- 
lar bar we use c/a = 0.1. For the dark matter bar 
c/a = 0.5. The other axial ratio b/a was varied for 
both ellipsoids to test its effect on the torque. We 
also varied the angle between the bars. The den- 
sity was chosen to fall exponentially with the dis- 
tance mimicking the decline of the density in bars 
measured in our N-body models. The ellipsoid 
representing the dark matter bar was twice shorter 
as compared with the stellar bar. We use Monte 
Carlo method to estimate the torque between the 
ellipsoids. We find that the torque r between the 
ellipsoids scales as r cx eoMCsteiiar sin(2A0), which 
is what one naively expects. Once the parameters 
of the ellipsoids are fixed, the model also provides 
absolute value of the torque. We used parame- 
ters of bars in the model D c b at 3 Gyr to roughly 
estimate the torque. Indeed, we find that the ex- 
pected torque between the bars and the measured 
torque in the model D c b are quite close (within 
30%). This is reassuring, but unfortunately, it is 
difficult to press this issue because matching the 
masses and the density distributions with bars in 
real simulations is difficult. Instead, we decided to 
use only the scaling relation for the torque. In this 
case we use the amplitude of the stellar bar A 2 as 
a measure of the bar strength. The amplitude of 
the expected torque is fixed to have the observed 
value at one moment in our models: at 3 Gyr. 
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We tested three of our models. In Figure 12 
we show the rate of angular momentum loss of the 
disk as a function of time (top panels) for mod- 
els D c b and Khb- In the bottom panels we show 
the evolution of the quantity r = A 2 eo sin(2A0) 
measured at r = Rd and normalized at 3 Gyr. A 2 , 
er>, and Acf> arc the amplitude of the disk bar, the 
cllipticity of the DM bar, and the phase between 
the bars, respectively, introduced in previous sec- 
tions. With the exception of the first ~ 1 Gyr 
(model D c b), the observed torque dL z /dt and the 
estimated torque between the bars are remarkably 
similar. After initial stage, at t 1 — 2.5 Gyr the 
torque declines. The second growth of the stel- 
lar bar at about 4 Gyr in model D c b produces an 
increase in dL z /dt at later epochs. 

We also find an agreement between the mea- 
sured and estimated torques for other models. We 
note that none of the quantities in the estimate 
of t can separately explain even qualitatively the 
time behavior of dL z /dt. Indeed, the decline in 
the torque from 1.5 Gyr to 3 Gyr is mostly due to 
decreasing of the angle between the bars. Yet, the 
steep increase in dL z /dt at 4 Gyr in model D c b 
is not related with the angle: the angle keeps de- 
creasing. The increase is mostly due to the growth 
of the amplitude of the stellar bar, which is accom- 
panied by an increase in the ellipticity of the dark 
mater bar (see Figure 8). In spite of general agree- 
ment of the torques, there are some deviations. At 
early stages (t < 1 Gyr) the scaled value of r is too 
small as compared with the measured torque in the 
cold model D c b- This is hardly surprising because 
the disk in this model develops instabilities very 
quickly. In addition, it was difficult to set the very 
central ~ 200 pc of the disk initially in equilib- 
rium. So, irregular non-axisymmetric waves are 
present almost from the very beginning. These 
waves result in ~ 1 — 2% loss of the disk angular 
momentum. At around 0.1 Gyr the disk develops 
strong spiral arms, which likely account for a large 
fraction of the torque at that time. As the spiral 
arms die and the bar grows, the torque starts to 
be dominated by the interaction between the dark 
matter and the stellar bars, and the r approxima- 
tion works better. 

Judging by the similarities between dL z /dt(t) 
and r(i), we find it plausible that the torque be- 
tween the stellar and the dark matter bars is re- 
sponsible for a large fraction of angular momen- 



tum transfer between the disk and the dark mat- 
ter. 
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Fig. 12. — The top panels show the specific torque 
dL z /dt experienced by the stellar component in 
models D c b and Khb with a strong and a moder- 
ate stellar bar, respectively The bottom panels 
show the evolution of an estimate of the torque 
due only to interaction between the stellar and 
the dark matter bars: t = eoA 2 sin(2A</>), where 
£d , A 2l and A(j> are the ellipticity of the halo bar, 
the amplitude of the stellar bar, and the lag angle, 
respectively. Torque dL z /dt is reasonably well re- 
produced by interaction between the bars during 
most of evolution. 

7. Discussion and conclusions 

Just as many other groups in the field (e.g., 
Athanassoula 2005a,b; Shen & Sellwood 2004; 
Debattista & Sellwood 1998) we find that our 
disk models develops strong, moderate, and weak 
bars. They show what seems to be a generic fea- 
ture of bar simulations: the growth in their am- 
plitude is accompanied with a decrease in their 
pattern speed. The role that the halo plays in this 
phenomenon appears to be crucial. The dark mat- 
ter ellipsoidal figure, the halo bar, that arises due 
to the presence of the disk and stellar bar, exerts a 
torque on the stellar bar that slows down this lat- 
ter. It is thus not surprising to find that the first 
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period of bar growth (the only one in some mod- 
els) coincides with a decline in the phase difference 
between the disk and the halo bar. Interestingly, 
an opposite trend seems to exist as well: bars ro- 
tate faster as they weaken. This is clearly seen in 
models K c b and Dh s - 

Disk particles located in the middle 2-8 kpc re- 
gion transfer angular momentum to the halo and 
to the outer disk region. In particular, this effect 
would explain the increase of the central stellar 
density seen, for example, in Figures 8-9 of Valcn- 
zuela & Klypin (2003) (see also Avila-Reese et al. 
2005). We have seen in section 4 that this disk 
density increase is accompanied by a correspond- 
ing central halo density increase. We apply the 
adiabatic contraction formalism to study whether 
the compression of the halo could be modeled by 
it. We find that the modified formalism by Gncdin 
et al. (2004) reproduces relatively well the inner 
density growth whereas the standard formula typ- 
ically overestimates the density increase. 

We compute the ellipticities of the halo bar 
in a plane parallel, ejj, and perpendicular, /, to 
the disk plane. The ellipticities decrease and the 
halo bar becomes more spherical as the radius in- 
creases. In some models, after an initial increase 
the flatness stays roughly constant or decreases 
during the disk evolution. On the other hand, in 
models with strong bars like model D c b, the el- 
lipticity €d has a second period of increase, which 
agrees with the second growth period of the stellar 
bar. As a result, at end of the evolution the in- 
ner dark matter halo ends up with a near prolate- 
like shape. Our results roughly agree with those 
mentioned in Athanassoula (2005a, b), these later 
being preliminary 

In all models, the inner halo (r SS 5kpc) is flat- 
tened, with the minor axis pointing perpendicular 
to the disk plane. The alignment of the minor axis 
with the disk axis is compatible with the results by 
Bailin et al. (2005). In cosmological simulations 
of disk galaxy formation Bailin et al. (2005) find 
that the the inner halo (r 0.1i? v ir) orients its mi- 
nor axis parallel to the disk axis, regardless of the 
orientation of the outer halo. As they explain, the 
misalignment between the two halo regions should 
be taken into account when modeling tidal streams 
in the halos of disk galaxies. 

Long et al. (1992) found that the presence 
of a stellar bar in the center of our Galaxy does 



not strongly affect the destruction rates of globular 
clusters. The results of Long et al. (1992) or of 
similar studies (e.g., Pichardo et al. 2004) might 
change in light of our findings: halo bars could 
enhance the dynamical effect of bars on globular 
clusters distribution. 

We use the same method of the tensor of inertia 
to study the evolution with time of the orientation 
of the halo bar relative to the disk bar. In all 
models, an alignment within ~ 10° is obtained 
after the initial bar growth. In some models, for 
example model Khb, the bars stay aligned during 
the whole evolution. It is interesting to see that 
in the model D c b the periods of alignment roughly 
agree with the two periods of stellar bar growth. 
In turn, the periods of bar growth are also the 
periods of higher angular momentum lost rate. 

As Valenzuela & Klypin (2003), we also find 
that the angular momentum transfer from the disk 
to the halo is rather modest: it is less than 10%. 
As it was shown by Valenzuela & Klypin, this is 
partly due to the fact that a significant fraction of 
the disk angular momentum is not lost but redis- 
tributed: outer regions absorb part of the angular 
momentum of intermediate regions. 
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